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Abstract 

IT) ' 

The solvability of The Dirac equation is studied for the exponential-type potentials with the 
pseudospin symmetry by using the parametric generalization of the Nikiforov-Uvarov method. 
The energy eigenvalue equation, and the corresponding Dirac spinors for Morse, Hulthen, and q 
deformed Rosen-Morse potentials are obtained within the framework of an approximation to the 
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spin-orbit coupling term, so the solutions are given for any value of the spin-orbit quantum number 
K = 0, or k 7^ 0. 
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I. INTRODUCTION 



The concept of pseudo-spin was constructed firstly in spherical nuclei [1, 2], and it is 
observed experimentally that the single particle levels labeled as pseudo-spin doublets are 
very close in energy [3] . The pseudo-spin doublets in nuclei are decomposed by using radial 
(n r ), orbital (£), and total angular momentum (j) quantum numbers as (n r ,£,j = £ + 1/2), 
and (n r — 1,£ + 2,j — £ + 3/2). The pseudo orbital angular momentum, £ = £ + 1, and 
the pseudo spin, s, quantum numbers give total angular momentum, j = £ + s, and the 
pseudo-spin doublets are degenerate with respect to pseudo spin, s [4]. The pseudo-spin 
doublets occur in nuclei, when the magnitude of scalar, V s (r), and vector, V v (r), potentials 
are nearly equal, with opposite sign, i.e., V s (r) ~ — V v (r) [5]. The pseudo-spin symmetry 
is studied based on Dirac equation in real nuclei, and shown that it is related with the 
competition between the centrifugal barrier, and pseudo-spin orbital potential [6]. The 
pseudo-spin concept is discussed in deformed nuclei [7], and exotic nuclei as well [3]. It is 
observed that the pseudo-spin symmetry is also an important one in the case of triaxiality 
[8]. 

The pseudo-spin orbital potential creates the splitting of the pseudo-spin doublets, and 
the pseudo-spin symmetry is exact symmetry in real nuclei, when the derivative of the 
difference between scalar, and vector potentials vanishes, but the above condition gives a 
good symmetry for exotic nuclei [3, 6]. The pseudo-spin symmetry is identified as a SU{2) 
symmetry of the Dirac Hamiltonian, under the condition that the sum of scalar, and vector 
potentials is equal to zero [9]. Recently, it is pointed out that the shape of the lower 
components of the Dirac spinor for the doublets is the same, when the pseudo-spin doublets 
are degenerate [4]. It is showed that it becomes possible to construct a map which relates 
the normal state (£, s) with the pseudo state (£, s) by applying of the helicity operator to 
the non-relativistic single-particle eigenfunction to understand of the mechanism to generate 
the pseudo-spin symmetry [10]. The pseudo-spin concept is discussed in the non-relativistic 
harmonic oscillator, and obtained that the condition between the coefficients of spin-orbit, 
and orbit-orbit terms in the case of non-relativistic single-particle Hamiltonian having the 
pseudo-spin symmetry is consistent with the result obtained relativistic mean-field theory 
[11, 12]. 

The pseudo-spin symmetry concept has found a great application area, especially in 
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nuclear theory. The identical bands observed in nuclei are explained by using the pseudo- 
spin symmetry [13]. The idea of the pseudo-spin has been considered to be useful to construct 
an effective shell- model coupling scheme [14]. Some features of nuclei, such as deformation, 
and super deformation can be explained in the context of the pseudo-spin symmetry [13, 15, 
16]. 

Recently, the solutions of the Klein-Gordon, and Dirac equations including the spin-orbit 
coupling term have been studied by many authors for different potentials, such as Morse 
potential [17-20], Poschl- Teller potential [21-23], Woods-Saxon potential [24], Eckart [25-28], 
harmonic oscillator [29, 30], three parameter diatomic molecular potential [31], and angle- 
dependent potential [32]. In the present work, we deal with the approximate solutions, 
and corresponding wave functions of the Dirac equation including spin-orbit coupling term 
under the exact pseudo-spin symmetry for Morse, Hulthen, and g-deformed Rosen-Morse 
potentials. We point out that the parametric generalization of the Nikiforov-Uvarov method 
can be applied to the Dirac equation with the above potentials, and the energy eigenvalue 
equation, and corresponding eigenfunctions can be obtained for the values of spin-orbit 
quantum number k = 0, or k ^ 0. 

II. THE DIRAC EQUATION WITH SPIN-ORBIT COUPLING 

The Dirac equation for a fermion with mass m moving in an external scalar, and vector 
potentials reads (h = c = 1) 

[a.p + P[m + V a (r)]+V v (r)]*(r)=E*(r), (1) 

where E is the energy of the particle, p is the three-momentum operator, and a, and [3 are 
the 4x4 Dirac matrices written in terms of2 x 2 Pauli matrices, and unit matrix. Under 
the consideration that the system has a spherical symmetry for which the potential fields 
depend on the radial coordinate, the quantum state of the particle is labeled by the quantum 
number set (n r , j,m, k), where m is the projection of the total angular momentum on the 
^-axis, and k = ±(j + 1/2) is the eigenvalues of the operator k = —(3(<7- L + 1) [33]. Here, 
n = —{j + 1/2) denotes the aligned spin (si/2,P3/2, etc.), and k = +(j + 1/2) denotes the 
unaligned spin (pi/2, ^3/2, etc.). The spherically symmetric Dirac wave function can than be 
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written in terms of upper, and lower components as [34] 



»(r)=i( /(r)[y " £l M (2) 

where f(r), and g(r) are the radial wave functions, Ye(9,(f>), and x are the spherical, and 
spin functions, respectively. Substituting Eq (2) into Eq. (1), we get the following radial 
Dirac equations 



(j r + £)/(r)-[M(r)+c]0(r) = O, (3) 
(J: - ^»(r)-[M(r)-c]/(r)=0. (4) 

where M(r) = m + V^(r), and e = E — V v (r). Using the expression for g(r) obtained from Eq. 
(3), and inserting it into Eq. (4), we have two second order differential equations including 
spin-orbit coupling term 



{ 55 - ^ - KM - ^}fM = {j^i^ - + 7>}/M • < 5 > 

{ £ - ^ - - e*]} 9 (r) = {-^ | K (r) + V«r)]( J; - =)} 9 (r) , (6) 

Under the condition of the exact pseudo-spin symmetry, i.e., ^[K(r) + V s (r)] = 0, or 
V v (r) + V s (r) = C = const., Eq. (6) gives 



{^2~ ^ +[m-E + C][V v (r)-V s (r)}}g(r) = [m 2 -E 2 + C(m + E)]g(r). (7) 

From the last equation, the energy eigenvalues depend also on the quantum number £ because 
of the relations given by k(k— 1) = £(£+!), and k(k+1) = £{£+!). So, the energy eigenstates 
with j — £± 1/2 are degenerate for £ ^ 0, which gives the situation of the exact pseudo-spin 
symmetry in the Dirac equation. 

In the present work, we intend to solve the last equation four different potentials, Morse, 
Hulthen, and g-deformed Rosen-Morse potentials, namely. The Dirac equation in Eq. (7) 
can not be solved exactly because of the spin-orbit coupling term. So, we use the Pekeris 
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approximation [35] to find the suitable expression instead of the spin-orbit coupling term in 
the case of Morse potential. This approximation makes possible to write the interaction term 
in the form of the Morse potential in terms of new parameters Di(i = 0, 1, 2). In the case of 
Hulthen, and g-deformed Rosen-Morse potential, we introduce the following approximation 
instead of the spin-orbit coupling term [22] 



1 a 2 e~ ar 
^ ~ (1 - e— ) 2 ' 

which gives a second order differential equation from Eq. (7) without (l/r 2 )-term. In all 
case of potentials, we use the parametric generalization of the Nikiforov-Uvarov method 
[36] , and so we prove that the parametric version of the method can be applied to the Dirac 
equation with Morse, Poschl- Teller, Hulthen, and g-deformed Rosen-Morse potential. 



III. NIKIFOROV-UVAROV METHOD 



By using an appropriate coordinate transformation, the Schrodinger equation is trans- 
formed into the following form 



where cr(s), a(s) are polynomials, at most second degree, and f (s) is a first degree poly- 
nomial. In the NU-method, the polynomial n(s), and the parameter k are required, and 
defined as 



n(s) = \ \a\s) - f(s)} ±J\ [a'(s) - f (*)]* ~ a(s) + ka(s), 



(10) 



and 



\ = k + 7r'(s), (11) 

where A is a constant. The function under the square root in the polynomial in n(s) in Eq. 
(10) must be square of a polynomial in order that ir(s) be a first degree polynomial, so the 
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derivative of n(s) is a constant, and this defines the constant k. Replacing k into Eq. (10), 
we define 



T (s) =f(s) + 27r(s). (12) 

where the derivative of r(s) should be negative [35], which let us know the choice of the 
solution. The hypergeometric type equation in Eq. (9) has a particular solution with degree 
n, if A in Eq. (11) satisfies 



X = X n = -„r>- W n - 1 ^\ „ = 0,1,2,... (13) 

To obtain the solution of Eq. (9) it is assumed that the solution is a product of two 
independent parts 



(14) 

where y(s) can be written as 



vM = ^ ) £ : W n (s)p{s)] 1 (is) 

where a n is normalization constant, and the function p(s) is the weight function, and should 
satisfy the condition 



^ P ( S ) + ,( S )^ = T(SMS ), (16) 
The other factor is defined as 



t)'(s) 7r(s) 



(17) 



<f)(s) a(s) 

In order to clarify the parametric generalization of the NU method, let us take the following 
equation, which represents a general form of the Schrodinger-like equation written for any 
potential, 
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d 2 ai - Ot 2 S d -£is 2 + £ 2 s ~ £3 
ds 2 s(l — 0:3s) ds [s(l — 0:3s)] 2 



When Eq. (18) is compared with Eq. (9), we get 



*(s) = 0. 



r(s) = «i - a 2 s ; cr(s) = s(l - a 3 s) ; ct(s) = -£is 2 + 6 s — 6 



3 • 



Substituting these into Eq. (10) 



(18) 



(19) 



7r(s) = «4 + 0:5s ± v/ («6 — fco^S 2 + («7 + fc)s + Ctg , 



(20) 



where the parameter set are 



(21) 



«4 = |(1 - on) , «5 = | («2 - 2a 3 ) , 
"6 = "5 + 6 > «7 = 2a 4 a 5 - £ 2 , 
«8 = «1 + £3 • 

In NU-method, the function under the square root in Eq. (20) must be the square of a 
polynomial [35]. This condition gives the roots of the parameter k, and they can be written 

as 



h,2 = —(017 + 2a 3 o; 8 ) ± 2 v /a^a^, 



(22) 



where the /c-values can be real or imaginary, and «g = 030:7 + a^as + olq . Different fc's lead 
to the different 7r(s)'s. For 



fc = -(a 7 + 2a 3 o; 8 ) - 2^a$a 9 , 



(23) 



7r(s) becomes 



7r(s) = a 4 + «5S - [(v 7 "^ + ttsv 7 ^) 5 ~ v 7 ^] > 



(24) 
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and also 



t(s) —ai + 2a 4 - (a 2 - 2a 5 )s - 2 [(y/ag + a 3y /c^)s - y/cts] ■ (25) 

Thus, we impose the following for satisfying the condition that the derivative of the function 
t(s) should be negative in the method 



r '( s ) = -(«2 - 2a 5 ) - 2( v /«^ + a 3v /o^) 

= -2a 3 - 2(^ + ^3^) <0- (26) 

From Eqs. (11), (12), (25), and (26), and equating Eq. (11) with the condition that A should 
satisfy given by Eq. (13), we find the eigenvalue equation 

a 2 n-(2n+l)a 5 + (2n + l)(y/ag~ + a 3 ^/ai) + n(n - l)a 3 

+ a 7 + 2a 3 a 8 + 2y/a 8 ag = 0. (27) 

We get from Eq. (16) 

p{s) = s ai °- 1 {l-a 3 s)^- ai0 - 1 , (28) 
and substituting into Eq. (15) gives 



y n (s)=pt°- 1 '^- ai0 - 1 \l-2a 3S ), (29) 

where 



«io =oii + 2a 4 + 2y/a$~ ; an = a 2 — 2a 5 + 2(y/a$~ + a 3y /a^) . (30) 
and p( a 'P\l — 2a 3 s) are the Jacobi polynomials. From Eq. (17), one gets 

<f>(s) = s^(l-a 3 sy ai *-^ , (31) 
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then the general solution \P(s) = <p(s)y(s) becomes 



V(s) = s ai *(l-a 3 s)- ai2 - ! 2?Pn 10 ^ ai ° 1) (l-2a 3 s) (32) 



where 



a 12 = a 4 + a/o^ ; a>i 3 — a 5 — (y/aQ + a 3v /a^) . (33) 



IV. BOUND-STATE SOLUTIONS 



A. Morse Potential 



We set the difference between scalar, and vector potentials in Eq. (6) as the Morse 
potential [37] given by 



V M orse(r) = D[e- 2ax - 2e- ax ] , (34) 

where D is the dissociation energy of the potential, a = ar , x = ^ — 1, r is the equilibrium 
distance, and the positive parameter a is the potential width. According to the Pekeris 
approximation, the spin-orbit coupling term can be written in terms of three new parameters 
D ,Di, and D 2 as 



K(K 2 1} ^ 7(A) + D 1 e~ ax + D 2 e~ 2ax ) , (35) 

where 7 = k(k — 1))tq (see Ref [17], for details). Substituting Eq. (34), and (35) into Eq. 
(6), we get 



^ + {ale' 2 - - ale— - a\)g{x) = , (36) 

where a\ = Tq^Dq + S m (m + E)], a\ = r^Di + 2DT, m ), and a| = r^DY^m — 7-D2) in the 
above equations. By using the new variable s = e~ ax , we obtain 
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?m + im + {al f-«3f-"jf ]9(s)=0 . (37) 

where ft 2 = 1/a 2 , and S m = m — E + C '. Comparing the Eq. (37) with Eq. (18), we obtain 
the following parameter set 



(X\ = 


1, 


-6 = 


« 2 = 


o, 


6 = ~a 2 2 (3 2 


a 3 = 


o, 


-6 = -a?/3 


(X4 = 


o, 


a 5 = 


a 6 = 


6, 


«7 = -6 


a 8 = 




"9 = 6 


a w = 


= 1 + 2^ 


an = 2a/^" 


Otyi = 


--Vis, 


«i3 = 



(38) 



Substituting the above expressions in Eq. (27), we obtain the energy eigenvalue equation of 
a fermion moving in the Morse potential under the exact pseudo-spin symmetry 

{2n+l)yJiD 2 -DZ m + /3r ( 7 £>i + 2L>E m ) 

+ 2Pr oy /[iD Q + S m (m + E)](^D 2 - DS m ) = , (39) 

and the lower spinor component from Eq. (32) 



X L 2 f^o + ^m +E) m .^ lD2 _ D ^ mS) m 

where L^\x) are the Laguerre polynomials. 

We give the upper spinor component from Eq. (4), by using Eq. (40) 
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+ L n -i (2^r \/7- D 2 - ^>S m s) 



x 



Pro^Do + S m (m + £) + — _ fag 



where we use some recursion relations of the Laguerre polynomials. 



B. Hulthen Potential 



The Hulthen potential reads 



(41) 



V(r) 



(42) 



where x — ^, and r denotes the spatial range [38, 39]. We choose the difference between 



scalar, and vector potentials as Hulthen potential potential, and using the approximation in 



Eq. (8), which becomes \ 



in the present case, we obtain from Eq. (6) 



d 2 g(r) j k(k — l)e 



dr 2 \r 2 (e x -l) 
Using the new variable s = e~ x , Eq. (43) becomes 



(43) 



d 2 g(s) 1 - s dg{s) 



ds 2 



s(l-s) ds [s(l-s)] : 



- rpl m (m + E) 



+ 
+ 



- k(k - 1) + 2r^S m (m + E) - r 2 12 m V 

- r 2 S m (m + E) + r 2 Z m Vo\ s 2 )g(s) = 



(44) 



Comparing with Eq. (18), we get 
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«3 



0:4 = , 

«6 = \ + 6 , 

aio = 1 + 2^ 

"12 = , 



-a = r 2 S m [-m - £ + Vo\ 

6 = r 2 S m [2(m + £) - Vb] - k(« - 1) 

= -r^S m (m + £) 
« 5 = -1/2 

«7 = -6 

"9 = 6-6 + 6 + 1 

an = 2 + 2(^1-6 + 6 + 1 + 

«i3 = -|-^/6-6 + 6 + i -v^s 



(45) 



The energy eigenvalue equation is written from Eq. (27) 



r ^/E m (m + E) [2n + 1 + ^Ak(k- 1) + 1 ] + (n + ^ )^/4k(k- 1) + 1 



+n(n + 1) + «(k - 1) + r 2 S m \/ + - = , 



(46) 



and corresponding lower, and upper Dirac spinors by using Eqs. (4) and (32), respectively, 



g ( s j = S r ^m{m+E) ^ _ s j i[i +v /4«(«-l)+l ] 
x ^,(2r ()v /S m (m+i?) iV /4k(/c-1)+1 )^ _ ^ 



f(s) = (r^y 1 s W Sm ( m+£ >(l - s )|[i+v /4 «(- 1 )+ 1 1 

x { [ - r oy /x m (m +~E) + -[1 + ^/4k(k - 1) + 1 ]^- 

« 1 ^,(2r 0v /S m (m+i;) , ^(/c-lHl ) ^ ^ , 

+ (n + 2r ^S m (m + E) + ^/4k(k - 1) + 1 + l)s 

(l+2r ov /S m (m+£) ,1 +V / K ( K -1)+1 ) 
X -r n _! (_! — AS) . 



(47) 



(48) 



where P^ ,m \x) are the Jacobi polynomials, and we have used some recursion relations of 
the Jacobi polynomials to obtain Eq. (48). 



-Deformed Rosen-Morse Potential 



The g-deformed Rosen-Morse Potential reads [40] 
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V(r) 



Vi 



,-2z 



- V 2 q- 



(49) 



1 + ge- 2 - " ^ (1 + ge- 2 -) 2 ' 
where x = ar, and g is the deformation parameter. 

We prefer to solve the Dirac equation for k = in Eq. (6) with the PT-symmetric version 
of the potential by setting to the difference between scalar, and vector potentials 



V 



PTi 



-2ix 



V 2 q- 



(50) 



1 + qe~ 2ix " ^ (1 + qe~ 2ix ) 2 ' 
We obtain the following equation under the above consideration by using the transformation 



s = — e 



-2ix 



d 2 g(s) 1-qs dg(s) 1 rr 2 r T/v v / . f m 

~d#- + 7(T^qT)— + ls(l-qs)] 2 { 6 ^m-Xmim + E)] 

+ 5 2 £ m [2g(m + E) - qV 1 + V 2 ]s - S 2 q 2 E m (m + E)s 2 }g(s) = , 
Comparing the Eq. (51) with Eq. (18), we get 



(51) 



«i = 1 , 
"2 = q, 
«3 = q, 

«4 = , 

"6 = ^+6, 

"8 = 6, 



-6=5VS m (m + P) 

G = 5 2 S m [2g(m + P)-g\/ 1 + \/ 2 ] 

-^ = 5 2 [VJl m -^ m (m + E)} 

a 5 = -q/2 

«7 = -6 

«9 = 6 " 96 + g 2 6 + 4 



"10 
«12 



1 + 2^, an = 2g + 2(J£i - ?6 + g 2 6 
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The energy eigenvalue equation is written from Eq. (27) 



(52) 



25y / S m (m + E - V,) [2n + 1 + ^ ^1 - (4V 2 X m )/q ] + (n + ^yl - (4V 2 E m )/ 9 

+n(n + l) + ^-5 2 £ m (\/ 1 + \/ 2 )/g = 0. (53) 

and corresponding lower, and upper Dirac spinors by using Eqs. (4) and (32), respectively, 
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x P fv^^.^(^) (1 _ 2gs) > (54) 
/(s) = 2iaE- 1 S V Sm ( fn+£ - Vl )(l- g s)3[ 1 +V 1 -( 4 ^ s ™)/9] 

x { [ - c^E^m + £ - + |[1 + ^/l - (4\/ 2 S m )/g ] T 



« ] (25 V /E m (m+E-Vl) )V /l-(4y 2 S m )/ ? ) 

+ g(n + 25y/X m (m + E-V 1 ) + ^1 - (4V 2 X m )/q + l) 

x ^VM^.iWHW/,)^ _ 2gs) (55) 



where <5 2 = l/(4a 2 ) in the above equations. 



V. CONCLUSION 



We have approximately solved the Dirac equation for the Morse, Hulthen, and g-deformed 
Rosen-Morse potentials with the exact pseudospin symmetry for arbitrary spin-orbit quan- 
tum number k. We have found the eigenvalue equation, and corresponding Dirac spinors 
in terms of Jacobi (or Laguerre) polynomials by using the parametric generalization of the 
NU-method within the framework of an approximation to the spin-orbit coupling term. We 
have observed that the parametric form of the NU method can be used to solve the Dirac 
equation with the above potentials. Our results for the cases of Morse potential is good 
agreement with the ones obtained in the literature. 
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